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Abstract 

A  model  and  an  analytical  solution  for  the  model  are  presented  for  the  resistance  of  the  polymer  electrolyte  membrane  of  a  H2/O2  fuel 
cell.  The  solution  includes  the  effect  of  the  humidity  of  the  inlet  gases  and  the  gas  pressure  at  the  anode  and  the  cathode  on  the  membrane 
resistance.  The  accuracy  of  the  solution  is  verified  by  comparison  with  experimental  data.  The  experiments  were  carried  out  with  a  Nafion  1 12 
membrane  in  a  homemade  fuel  cell  test  station.  The  membrane  resistances  predicted  by  the  model  agree  well  with  those  obtained  during  the 
experiments. 
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1.  Introduction 

The  pressing  need  for  a  clean,  efficient  and  reliable  power 
source  has  led  to  the  rapid  development  of  fuel  cell  technology. 
Fuel  cells  employing  a  polymer  electrolyte  membrane  (PEM)  are 
widely  regarded  as  a  potential  power  source  for  many  automo¬ 
tive  and  stationary  applications  [1,2].  These  polymer  electrolyte 
membrane  fuel  cells  (PEMFCs)  consist  of  the  membrane  elec¬ 
trode  assembly  (MEA)  and  gas  channels.  The  MEA  is  made 
up  of  two  gas  diffusion  layers  (GDLs),  two  carbon  supported 
catalyst  layers  and  PEM  [3].  The  GDLs  allow  the  simultaneous 
transport  of  the  gas  (H2  or  O2)  and  water  from  the  gas  channels 
to  the  catalyst  layers  and  also  serve  as  the  current  collectors. 
The  catalyst  layers  (commonly  Pt  or  the  Pt  alloys)  provide  the 
active  sites  required  for  the  electrochemical  reactions.  During 
the  MEA  preparation  the  catalyst  layer  is  either  applied  on  to 
the  GDL  or  on  to  the  membrane,  followed  by  the  addition  of  the 
other  component  [4].  The  core  of  the  PEMFC  is  the  membrane 
that  is  sandwiched  between  the  catalyst  layers.  The  PEM  con¬ 
ducts  protons  from  the  anode  to  the  cathode  and  separates  the 
anode  and  the  cathode  gases  from  each  other.  A  detailed  review 
of  some  of  these  proton  conductors  and  their  application  to  fuel 
cells  is  presented  by  Kreuer  [5],  The  most  commonly  used  mem¬ 
brane  in  the  PEMFC  is  Du  Pont’s  Nafion  [6].  Its  properties  and 
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applications  to  fuel  cells  are  extensively  reported  in  the  literature 
[7-10], 

The  conductivity  of  the  PEM  is  a  key  factor  for  the  perfor¬ 
mance  of  the  fuel  cell  system.  There  are  microscopic  models  in 
the  literature  that  discuss  the  membrane  conductivity  in  detail 
[11-18].  The  conductivity,  to  a  large  extent  depends  on  the  water 
content  of  the  membrane  [19].  The  macroscopic  models  avail¬ 
able  in  the  literature  are  either  the  diffusive  flow  type  models 
[20-22]  or  the  hydraulic  flow  type  models  [23,24].  The  diffu¬ 
sion  based  models  employ  an  empirical  expression  developed  by 
Springer  et  al.  [20]  or  its  modified  form  for  the  water  content  in 
the  membrane.  The  model  developed  by  Bernardi  and  Verbrugge 
[23]  forms  the  basis  for  most  of  the  hydraulic  flow  type  mod¬ 
els.  Apart  from  these  models,  there  are  also  other  macroscopic 
models  that  employ  chemical  potential  as  the  driving  force 
instead  of  specifying  the  flow  type  [25,26].  Recently  a  two-phase 
model,  which  includes  the  simultaneous  presence  of  liquid  and 
vapor  in  the  membrane,  was  developed  by  Weber  and  Newman 
[27]. 

The  present  study  is  based  on  the  hydraulic  flow  type  model 
developed  by  Eikerling  et  al.  [28],  which  includes  a  relation¬ 
ship  between  water  content  and  the  microscopic  structure  of 
the  membrane.  The  objective  of  this  study  is  to  predict  the  ionic 
resistance  of  the  membrane  using  a  simple  model  and  to  validate 
it  using  experimental  data.  An  attempt  is  made  to  introduce  some 
simplifying  assumptions  while  trying  to  maintain  the  physical 
insight  provided  by  the  original  model. 
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Nomenclature 

Cw 

molar  concentration  of  water  (mol  m-3) 

f 

fraction  of  hydrophobic  pores  in  the  GDL 

F 

i 

Faraday’s  constant  (C  mol-1) 
current  density  (A  cm-2) 

I 

current  (A) 

Km 

permeability  of  the  PEM  (m2) 

< bs 

absolute  permeability  of  the  GDL  (m2) 

is  m 
^-abs 

absolute  permeability  of  the  PEM  (m2) 

<  i 

relative  permeability  of  the  GDL 

Ld 

length  of  the  GDL  (m) 

Lm 

length  of  the  PEM  (m) 

< 

inlet  molar  flux  of  hydrogen  (mol  m-2  s) 

< 

inlet  molar  flux  of  oxygen  (mol  m-2  s) 

..darcy, 7 

JVw 

molar  flux  of  water  due  to  the  Darcy’s  flow  in  the 
region  j  (mol  m-2  s) 

NdrdS-j 

tvw 

molar  flux  of  water  due  to  the  electro-osmotic  drag 
in  the  region  j  (mol  m-2  s) 

A/in’a 

tvw 

molar  flux  of  water  at  the  anode  GDL  inlet 
(mol  m-2  s) 

Afin.c 

JVW 

molar  flux  of  water  at  the  cathode  GDL  inlet 
(mol  m-2  s) 

Nm 

tvw 

molar  flux  of  water  in  the  PEM  (mol  m-2  s) 

<°d 

molar  flux  of  water  produced  (mol  m-2  s) 

Pc 

capillary  pressure  in  the  region  j  (N  m-2) 

Pg 

gas  phase  pressure  in  the  region  j  (N  m-2) 

Pt iq 

liquid  phase  pressure  in  the  region  j  (N  m-2) 

Pw 

partial  pressure  of  water  (N  m-2) 

J 

rc 

capillary  radius  of  the  region  j  (m) 

rmd 

mean  pore  size  of  the  GDL  ( m ) 

R 

resistance  of  the  membrane  (mQ  cm-2) 

R* 

resistance  of  the  membrane  at  the  reference  con¬ 
dition  (mfl  cm-2) 

*w 

water  content  in  the  anode  GDL  (m3  m-3) 

Jw 

water  content  in  the  cathode  GDL  (m3  m-3) 

T 

absolute  temperature  (K) 

V(r) 

integral  pore  size  distribution  function  for  the 
membrane  (m3  m-3) 

wd 

mole  fraction  of  water  leaving  the  anode  humidi¬ 
fier  along  with  hydrogen  gas 

wc 

mole  fraction  of  water  leaving  the  cathode  humid¬ 
ifier  along  with  oxygen  gas 

X 

geometrical  coordinates  for  membrane  (m) 

f 

geometrical  coordinates  for  anode  GDL  (m) 

f 

geometrical  coordinates  for  cathode  GDL  (m) 

Greek  symbols 

a(r) 

pore  size  distribution  function  for  the  GDL 

y 

surface  tension  of  water  (N  m-1) 

s 

porosity  of  the  GDL 

p 

potential  drop  across  the  membrane  (V) 

9> 

contact  angle  of  the  region  j  (radian) 

k  specific  conductivity  of  the  membrane 

(1  f2  cm-1) 

k  Specific  conductivity  of  the  membrane  at  the  ref¬ 

erence  condition  (1  Q  cm- 1 ) 

X  water  content  of  the  membrane  ( m3 m-3) 

Xs  water  content  of  the  saturated  membrane 

(m3  m-3) 

/i  viscosity  of  water  (Nm-2s) 

$  electro  osmotic  drag  coefficient  of  the  membrane 

a  standard  deviation  of  the  pore  size  distribution  in 

the  GDL  (m) 

(p i  effective  porosity  of  the  hydrophilic  pores  in  the 
GDL 

(p2  effective  porosity  of  the  hydrophobic  pores  in  the 
GDL 

Subscripts 
abs  absolute 

atm  atmospheric 

C  capillary 

g  gas 

liq  liquid 

rel  relative 

S  saturation 

W  water 

Superscripts 
a  anode 

c  cathode 

d  GDL 

drag  electro-osmotic  drag 

darcy  Darcy’s  flow 

in  inlet 

m  membrane 


2.  Experiment 

2.1.  Membrane  electrode  assembly  (MEA)  preparation 

Pt  catalyst  ink  with  75wt.%  catalyst  and  25wt.%  Nafion® 
(dry  solids  content)  was  prepared  with  commercially  available 
40.2  wt.%  Pt  on  Vulcan  XC-72  catalyst  (E-TEK,  De  Nora  North 
America  Inc.,  NJ)  and  perfluorosulfonic  acid-PTFE  copolymer 
(5%,  w/v,  Alfa  Aesar,  MA).  Isopropyl  alcohol  (99%,  v/v,  VWR 
Scientific  Products)  was  used  as  the  thinning  solvent.  The  ink 
was  mixed  well  for  at  least  8  h.  ELAT  electrodes  (E-TEK,  De 
Nora  North  America  Inc.,  NJ)  coated  with  carbon  were  used  as 
gas  diffusion  layers  (GDLs).  The  GDLs  were  cut  into  10  cm2 
pieces.  The  catalyst  ink  was  then  sprayed  onto  the  GDLs,  air 
dried  for  30  min  and  then  dried  at  1 10  °C  for  10  min  to  evapo¬ 
rate  any  remaining  solvent.  The  process  was  repeated  until  the 
targeted  loading  was  achieved.  Both  the  anode  and  the  cathode 
sides  had  a  loading  of  0.5  mg  cm-2  of  Pt.  The  catalyzed  GDLs 
were  then  bonded  to  a  pretreated  Nafion®  1 12  membrane  by  hot 
pressing  at  140  °C  for  2  min  at  500  psig  to  make  the  MEAs.  The 
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MEA  was  assembled  into  the  fuel  cell  with  single  channel  ser¬ 
pentine  flow  field  plates  obtained  from  Fuel  Cell  Technologies 
in  which  4.25  mm  shims  were  added  to  the  assembly  to  ensure 
a  gap  thickness  that  would  provide  the  optimum  compression 
of  the  GDL.  The  assembly  was  compressed  at  200-325  psia 
(1.37-2  x  106Nm-2)  for  5  min.  The  cell  was  tightened  with 
eight  bolts  to  target  a  torque  of  20  in.-lbs  (2.25968  Nm)  on  each. 
The  cell  was  allowed  to  rest  for  25  min,  during  which  time  the 
necessary  torque  and  compression  corrections  were  applied  so 
the  final  pressure  of  200-325  psia  (1.37-2  x  106Nm-2)  and 
bolt  torque  of  20  in.-lbs  (2.25968  Nm)  were  achieved. 

2.2.  Measurements 

High  purity  H2  and  ultra  high  purity  O2  were  used  as  the 
fuel  and  oxidant,  respectively.  The  flow  rates  were  maintained 
at  200%  of  the  stoichiometric  requirements.  The  reactant  gases 
were  fully  humidified  (>100%  RH)  unless  otherwise  speci¬ 
fied.  The  cell  temperature  was  maintained  at  70  °C  throughout 
the  experiments,  while  the  anode  and  the  cathode  humidifiers 
were  maintained  at  75  and  70  °C,  respectively.  The  gas  flow 
rates  were  regulated  with  mass  flow  controllers.  The  fuel  cell 
was  tested  with  a  test  station  made  at  University  of  South  Car¬ 
olina.  The  flow  diagram  of  the  test  station  is  shown  in  Fig.  1. 
Before  collecting  the  performance  data,  the  cell  was  operated 
at  0.4  V  at  the  required  stoichiometric  ratios  of  the  gases,  for  a 
minimum  of  24  h.  The  polarization  data  were  collected  in  the 
constant  voltage  mode.  The  data  were  collected  with  a  volt¬ 
age  increment  of  30  mV.  An  open  circuit  rest  time  of  2  min 
was  allowed  before  collection  of  each  data  point.  The  contact 
resistances  were  made  negligible  by  compressing  the  MEA  to 
achieve  the  required  pinch  (0.33  mm)  specified  to  neglect  contact 
resistances  [29].  Hence  the  high  frequency  (1kHz)  resistance 


measured  (d)?/d/  =  d(/R)/d/)  using  a  Hewlett-Packard  milliohm 
meter  (Model  4328A)  is  mainly  attributed  to  the  membrane 
resistance  alone.  The  polarization  curves  and  the  resistance  mea¬ 
surements  were  recorded  using  a  personal  computer  interfaced 
with  LAB  VIEW  (National  Instruments).  The  electronic  load  for 
the  fuel  cell  system  was  Agilent  6060B  which  can  handle  a  cur¬ 
rent  up  to  60  A. 

3.  Model  description 

In  this  section  the  governing  equations  for  the  PEM  are  sim¬ 
plified  based  on  the  assumptions  stated  below.  This  provides 
an  analytical  solution  for  the  ionic  resistance  of  the  membrane. 
The  resistance  obtained  using  this  expression  is  compared  with 
experimental  data.  Further  the  analytical  solution  is  validated  by 
comparing  the  results  with  the  numerical  solution  developed  in 
a  subsequent  section. 

3.1.  Fuel  cell  system  and  basic  assumptions 

The  fuel  cell  model  considered  here  is  a  one-dimensional 
system  consisting  of  three  regions.  The  system  contains  an  anode 
GDL  and  a  cathode  GDL  in  addition  to  the  PEM. 

The  anode  and  the  cathode  catalyst  layers  are  not  considered 
explicitly  and  they  are  considered  to  be  very  thin  layers  of  neg¬ 
ligible  thickness  present  at  the  respective  GDL/PEM  interfaces. 
The  electrochemical  reactions  and  the  phase  changes  are  consid¬ 
ered  to  occur  at  these  GDL/membrane  interfaces.  The  schematic 
diagram  of  the  system  is  shown  in  Fig.  2.  The  basic  assumptions 
used  in  the  model  are  given  below: 

(i)  The  system  is  considered  to  be  isothermal. 

(ii)  The  system  is  considered  to  be  at  steady  state. 


Mass  flow  meter 


Mass  flow  meter 


Fig.  1.  Schematic  representation  of  the  fuel  cell  experimental  setup. 
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Fig.  2.  Schematic  representation  of  the  fuel  cell  model  system. 
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(iii)  Both  the  anode  GDL  and  the  cathode  GDL  are  considered 
to  be  identical  porous  media. 

(iv)  The  PEM  in  this  case  is  Nafion,  a  perfluorosulfonic  acid 
membrane  with  a  PTFE  backbone  [10].  The  PEM  is  con¬ 
sidered  to  have  a  porous  structure.  The  sulfonic  acid  group 
separates  itself  from  the  polymer  matrix  and  forms  clusters 
when  hydrated  with  water.  These  ionic  clusters  grow  and 
form  interconnecting  channels,  which  expand  to  accom¬ 
modate  the  liquid  water  that  is  in  equilibrium  with  the 
membrane  [13].  This  gives  rise  to  porous  structure. 

(v)  The  pores  in  all  the  three  regions  of  the  fuel  cell  system  are 
treated  as  cylindrical  capillaries. 

(vi)  The  capillary  framework  is  used  for  water  transport  in  all 
the  regions. 

(vii)  The  water  vapor  brought  along  with  the  anode  and  the 
cathode  gas  is  assumed  to  condense  into  liquid  water  at 
the  respective  GDL/membrane  interface.  There  is  no  phase 
change  considered  in  the  GDL  regions. 

(viii)  The  water  removal  from  the  system  is  considered  to  be  in 
the  liquid  phase. 

(ix)  The  membrane  is  assumed  to  be  in  equilibrium  with  only 
liquid  water.  No  gas  cross  over  is  considered  across  the 
membrane. 


There  is  no  net  water  flux  across  the  membrane  (see  Appendix 
A). 

These  assumptions  allow  us  to  treat  liquid  water  transport  in 
each  region  separately  using  capillary  phenomena  [28].  In  the 
context  of  capillary  framework,  the  characteristic  pore  size  is 
taken  as  the  radius  of  the  capillary  meniscus.  The  capillary  radius 
( rc )  is  related  to  the  capillary  pressure  (pc)  and  the  wetting 
properties  of  the  medium  are  given  by  the  Laplace’s  expression 
[30]: 


i  2ycos(@'7) 
rc  = - t - ,  j  —  a,  c,  m 

Pc 


(1) 


where  y  is  the  surface  tension  of  water  and  (P  is  the  contact 
angle  between  water  and  the  medium  j.  The  capillary  pressure 
(pg)  is  the  result  of  two  opposing  forces  acting  in  the  membrane, 
namely  the  gas  phase  pressure  (pJg)  and  the  liquid  water  pressure 

(Pi iq); 

p'c  =  Pg  ~  Ph q,  j  =  a,  c,  m  (2) 

In  the  following  section  the  model  of  the  membrane  system  is 
described  in  detail  and  an  analytical  solution  for  the  membrane 
ionic  resistance  is  obtained. 

3.2.  Membrane 

3.2.1.  Governing  equations 

The  anode/membrane  interface  is  chosen  as  the  origin  of  the 
coordinate  system  for  the  membrane  model  (see  Fig.  2).  The 
transport  of  the  protons  across  the  membrane  takes  place  in 
the  presence  of  water.  The  proton  drag  across  the  membrane 
is  described  by  Ohm’s  law  [31]: 

i  =  (3) 

ax 

rLm  dx 

V  =  -i  77 — ^  (4) 

Jo  K(X(pc)) 

where  i]  is  the  potential  drop  across  the  membrane,  K().(pc))  the 
specific  conductivity,  which  is  the  function  of  the  water  content, 
X(pc )  [28],  which  in  turn  is  a  function  of  the  capillary  pressure, 
pc-  A  linear  dependence  of  conductivity  on  water  content  is 
assumed  [32]: 

k(Hpc))  =  k*Kpc)  (5) 

where  k"  is  the  specific  conductivity  of  the  membrane  at  the 
reference  condition.  In  this  study  the  open  circuit  operating  con¬ 
dition  is  chosen  to  be  the  reference  condition. 
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The  flux  of  water  in  the  membrane  is  comprised  of  two  com¬ 
ponents.  One  is  due  to  the  electro-osmotic  drag,  i.e.,  the  water 
dragged  in  along  with  the  proton.  This  part  of  the  water  flux, 


Nfdg  is  given  by  [27]: 

^drag  _ 

"w  -  F 


(6) 


Here  i/F  is  the  proton  flux  across  the  membrane  and  q  is  the 
electro-osmotic  drag  coefficient,  which  is  the  number  of  moles 
of  water  dragged  along  with  each  mole  of  the  proton.  The  other 
part  is  due  to  the  Darcy’s  flow  [3]: 


darcy  _ 

ivw  — 


cwKm(X(pc))  dpiiq 
fi  dx 


(7) 


Hence  the  total  water  flux  inside  the  membrane  is  given  by: 


Arm  _  «£  cwKm(X(pc ))  d pnq 
Nw~F  dT 


(8) 


The  electro-osmotic  drag  coefficient  (£)  is  assumed  to  be  a  con¬ 
stant  [27].  For  the  permeability  of  the  membrane,  K'"(X(pc)), 
we  adopt  a  linear  dependence  of  water  content,  which  is  a  rea¬ 
sonable  approximation  for  thin  membranes  (50  p,m): 


Km(X(Pc ))  =  K™hsX(pc) 


(9) 


where  A'fT  is  the  absolute  permeability  of  water  in  the  mem¬ 
brane. 

Using  (2)  and  (8)  Eq.  (7)  can  be  written  as 


Nm  — 

7Vw  — 


«£  cwKtf>sWPc)  ( d(-pc)  .  dpg 


P 


dx 


dx 


(10) 


Since  we  have  assumed  the  gas  crossover  to  be  negligible, 
the  gas  pressure  profile  across  the  membrane  is 


dPg 

dx 


Pa~Pc 


APg 

Lm 


(11) 


Hence  Eq.  (10)  yields  the  net  water  flux  across  the  membrane, 
which  is  written  as 


Nm  — 

iVw  — 


Cw^s7.(pc)  ( d(-pc)  .  APg 


P 


dx 


(12) 


Since  the  net  water  flux  across  the  membrane  is  assumed  to 
be  zero  for  the  all  the  operating  conditions  considered  in  this 
study  (see  Appendix  A),  Eq.  (12)  becomes: 


d(-pc)  _  i%P _ Apg 

dx  “  cwK™hsX(pc)F  Lm 


(13) 


3.2.2.  Membrane  resistance 

Eqs.  (4)  and  (13)  form  the  set  of  governing  equations  for  the 
transport  of  species  (proton  and  water,  respectively)  in  the  mem¬ 
brane.  These  two  equations  are  combined  to  yield  the  following 
expression  for  potential  drop  across  the  membrane: 


rp  C2 


n  = 


d(-pc) 


variable  of  integration  from  pc  to  rc  and  utilizing  Eq.  (1),  Eq. 
(14)  becomes: 


i]  =  i2y  cos(6m)R* 
where 

o  (/£/  F)pLm 


d(r c) 


Irci  rc(P  -  APg^rc)) 


cw  K™ 


abs 


and 


R*  = 


(15) 


(16) 


(17) 


Comparing  Eq.  (15)  with  Ohm’s  law  i.e.,  Eq.  (4),  we  have  an 
expression  for  the  membrane  resistance: 


V 


rr c2 


d(rc) 


(18) 


=  —4  =  — 2y cos(0  )/?  /  -T- 

i  J  rc  1  ~  ApgX(rc)) 

In  Eq.  (18)  the  membrane  resistance  R  is  a  function  of  the 
current  density  ( i )  and  water  content  in  the  membrane  (/,(/'(-)). 
In  order  to  evaluate  the  membrane  resistance,  the  water  content 
of  the  membrane  (/,)  as  a  function  of  capillary  radius  (rc)  is 
necessary. 

3.2.3.  Expression  for  water  content  of  the  membrane 

The  expression  for  water  content  in  the  membrane  is  obtained 
by  virtue  of  the  formation  of  the  porous  structure  and  the  defini¬ 
tion  of  water  content  in  the  membrane.  As  mentioned  before 
the  sulfonic  acid  group  present  in  the  membrane  side  chain 
forms  clusters  and  interconnecting  channels  in  the  presence  of 
water.  When  liquid  water  enters  the  membrane,  these  channels 
interconnecting  the  clusters  expand.  The  water  content  in  the 
membrane  is  the  fraction  of  the  channels  that  has  expanded  to 
accommodate  the  liquid  water  [28].  It  is  determined  by  the  inte¬ 
gral  pore  size  distribution  function,  V(rc),  normalized  to  the 
saturation  water  content  for  the  membrane  (/,s),  which  is  equiv¬ 
alent  to  the  maximum  volume  of  the  pores  in  the  membrane 
[33]: 


Hrc)  = 


V(rc) 
^ s 


(19) 


In  the  case  of  Nation  112  the  maximum  water  content  is 
0.44  cm3  cm-3  [33]. 

3.2.4.  Analytical  solution  for  the  membrane  resistance 

A  straight  line  is  fitted  to  the  integral  pore  size  distribution 
data  [33]  ( V(/‘c)  versus  rc)  of  Nation  1 12.  Using  this  relation 
and  Eq.  (19)  we  get  the  relation  between  the  water  content  in 
the  membrane  and  the  pore  radius  as 


X(rc)  =  8607. 3rc  +  0.6 


(20) 


/m  (K*X(pc))m/F)p))/(cwKfhsX(pc))  -  (APg)/(Z™)} 


(14) 


where  pci  and  pc2  are  the  capillary  pressures  of  the  anode  and 
the  cathode  side  of  the  membrane,  respectively.  Changing  the 
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Using  (18)  and  (20),  we  arrive  at  an  analytical  expression  for 
the  membrane  resistance: 


R  =  —2ycos(0m)R* 


In 


(rC2-rci)8607-3A^2 


[0  -  Apg(8607.3(rC2  -  rCi)  +  0.6)]8607  3A^2 


1 

S(rC2  -  rCi) 


(21) 


where 

8  =  p-  0.6Apg  (22) 

In  Eq.  (2 1 ) ,  rc  i  and  rc2  are  the  capillary  radii  at  the  membrane 
GDL  boundaries,  which  are  estimated  using  the  GDL  model 
(developed  in  a  subsequent  section). 


The  water  content  of  the  GDL,  is  related  to  the  capillary 
pressure,  pJc  and  the  properties  of  the  GDL  by  the  following 
expression  (see  Appendices  A  and  B): 


^1  +  (1  -  2<fe)erf 


3.2.5.  Numerical  solution  for  the  membrane  resistance 

To  estimate  the  accuracy  of  the  analytical  expression  devel¬ 
oped  in  the  previous  section,  we  have  also  evaluated  the  mem¬ 
brane  resistance  using  a  numerical  model.  In  order  to  accomplish 
this  task  the  actual  integral  pore  size  distribution  data  for  Nation 
112  [33]  is  correlated  with  the  capillary  radius,  rc,  and  using 
Eq.  (19)  we  get  the  membrane  water  content  as  the  function  of 
the  capillary  radius: 


k(rc)  =  0.978 


8.04  x  10 


-9 


3.09  x  10 


-17 


4.57  x  10 


rc 

-26 


1.15  x  10“  35  1.34  x  10 


-44 


(23) 


Eq.  (23)  is  used  along  with  Eq.  ( 1 8)  and  the  GDL  model  to  obtain 
a  numerical  solution  for  the  PEM  resistance. 


3.3.  GDL  model 

The  GDL  is  considered  to  be  a  porous  medium  comprising 
of  a  mixture  of  hydrophilic  and  hydrophobic  pores  [34] .  The 
liquid  water  transport  inside  the  GDL  is  considered  to  be  only 
by  capillary  phenomena.  The  gas  phase  pressure  is  considered 
to  be  uniform  in  each  GDL  region.  Due  to  the  presence  of  gas 
in  the  medium,  the  flow  of  water  is  similar  to  that  of  the  flow  in 
an  unsaturated  medium  [  34] . 

The  water  flux  flowing  away  from  the  membrane  in  the  GDL 
is  described  by  Darcy’s  flow  for  unsaturated  medium  [34] ,  which 
is  given  as 


'  ln(— ((2y  cos  6d)/p())  -  ln(rmd)' 

ct-v/2 


J  =  a,c 
(27) 


Here  r/j?  is  the  porosity  of  the  hydrophobic  pores  in  the  GDL,  y 
the  surface  tension  of  water,  9d  the  contact  angle  made  by  water 
with  the  GDL,  rmd  the  mean  pore  size  in  the  GDL  and  o  is  the 
standard  deviation  of  the  pore  size  distribution  function  of  the 
GDL  material. 

The  liquid  pressure  or  capillary  pressure  profile  in  the  anode 
and  cathode  GDL  is  solved  using  a  simple  water  balance  in  the 
respective  region.  Since  the  water  flux  is  assumed  to  be  uniform 
in  the  GDL  and  no  phase  change  is  considered,  the  inlet  flux  of 
water  in  each  region  is  the  same  as  the  flux  of  water  described  by 
Darcy’s  flow  (away  from  the  membrane  towards  the  respective 
gas  channels)  for  that  region.  In  the  cathode  region  (cathode 
GDL/membrane  boundary)  there  is  additional  water  produced 
due  to  the  electrochemical  reduction  of  oxygen.  Hence  we  have: 

for  the  anode  region  :  A^.dCy'a  =  A^'d  (28) 

for  the  cathode  region  :  A^J.dcy'c  =  +  N^od  (29) 

The  feed  coming  out  of  the  anode  and  the  cathode  humidifiers 
is  considered  to  be  a  binary  mixture  of  the  water  vapor  and 
H2/O2  gas,  respectively.  Fick’s  law  of  diffusion  [35]  is  used  to 
calculate  the  flux  of  water  leaving  from  the  humidifier  along 
with  the  gases: 

for  the  anode  :  <’a  =  \  (30) 


A,dracy,;  _ 
lyW  ~  ~ 


CW^abs^el  d^q 


J  =  a,c 


(24) 


p  dy7 

The  relative  permeability  KJ.  is  assumed  to  be  a  cubic  func 

™  c 

w 


tion  of  water  content  in  the  medium  S^,  as  shown  by  [34],  Hence 


we  have: 


J 


Afdracy J  _  cW^abs(^7)  ^Fliq 

"  p  d^ 

Combining  Eqs.  (2)  and  (25)  we  have: 


a,  c 


n: 


dracy.j  _  cW^bs(^)  dPC 

p  d  yJ  ’ 


W 


]  —  a,c 


(25) 


(26) 


According  to  Faraday’s  law  [31],  N^2  —  1/2 F  is  the  the¬ 
oretical  molar  flux  requirement  of  hydrogen.  Hence  Eq.  (30) 
becomes: 


ATin,a  _  1  l  w 

tvw  —  — 


IF  V  1  -  ws 


and  for  the  cathode  :  = 


-N, 


02 


(31) 

(32) 


where  Ng2  =  i/4Fis  the  theoretical  molar  flux  requirement  of 
oxygen.  Hence  Eq.  (32)  becomes: 


Nm’c  -  - 


4F(l  -  wc) 


(33) 
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Along  with  this  the  flux  for  water  produced  (according  to 
Faraday’s  law  [31])  in  the  cathode  which  is  written  as 


NPmd  - _ 

Nw  ~  2  F 

Thus  water  flux  in  cathode  GDL  region  becomes: 


(34) 


Wln'c  -  _ _ 

iVw  — 


2F  4F(l  -  wc) 


(35) 


Using  Eqs.  (26),  (28)  and  (31),  the  governing  equation  for 
the  anode  region  becomes: 


dpi  =  (i/2F)(wa/(l  -  wa))p. 

^  cwKadbs(^) 


(36) 


Using  Eqs.  (26),  (29)  and  (35),  the  governing  equation  for 
the  cathode  region  becomes: 


dPCC 

dvc 


((//2 F)  +  iwc / (4 F([  -  wc)))ii 


(37) 


At  both  the  anode  and  the  cathode  GDL  inlets  (GDL/gas  channel 
boundary),  the  gas  and  the  liquid  pressure  continuity  is  assumed. 
Thus  we  have: 


Pg  =  Pdg  (38) 

and 

A  =  A j  (39) 

Since  there  is  no  liquid  in  the  gas  channel  we  let  the  liquid 
pressure  as  zero.  The  capillary  pressure  at  the  GDL/gas  channel 
boundary  (GDL  inlet)  is  given  by 


dan  _  „d,in 

PC  ~  P g 


d,in 

Fliq 


Pc 


Using  (38)  and  (39),  Eq.  (40)  becomes: 
Pgn-° 


(40) 

(41) 


This  is  used  as  the  boundary  conditions  to  solve  Eqs.  (36)  and 
(37)  to  obtain  the  pressure  profiles  at  the  anode  and  the  cath¬ 
ode  regions,  respectively.  These  differential  equations  along 
with  their  respective  boundary  conditions  are  solved  separately 
(dsolve,  numeric)  using  MAPLE  software.  The  capillary  radius 
at  a  membrane/GDL  boundary  is  calculated  using  the  capillary 
pressure  obtained  at  that  boundary,  and  Eq.  (1).  These  capillary 
radii,  i  and  rc2  obtained  at  the  membrane  boundaries  are  used 
in  the  calculation  of  the  membrane  resistance  in  Eq.  (20). 


Fig.  3.  Polarization  curves  under  different  operating  conditions,  (a)  Fully  humid¬ 
ified  gases  with  no  gas  pressure  at  the  anode  and  the  cathode,  (b)  Fully  humidified 
gases  with  cathode  gas  pressure  of  2  atm.  (c)  Fully  humidified  gases  with  the 
anode  gas  pressure  of  2  atm.  (d)  Dry  cathode  gas  (O2)  with  no  gas  pressure  at 
the  anode  and  the  cathode,  (e)  Dry  anode  gas  (H2)  with  no  gas  pressure  at  the 
anode  and  the  cathode. 


almost  to  the  same  degree  as  that  of  the  anode  at  lower  current 
density  region  (up  to  0.8  A  cm-2),  while  beyond  that  there  is 
a  gradual  increase  in  the  performance.  When  dry  oxygen  gas 
is  used  as  the  oxidant  the  performance  decreases.  There  is  a 
dramatic  decrease  in  the  performance  of  the  fuel  cell  system 
when  dry  FL  gas  is  used  as  the  anode  feed.  The  loss  in  voltage 
drop  across  the  membrane  is  nearly  100  mV. 


4.2.  Membrane  resistance 


4.2.1.  Resistance  measured  experimentally 

When  the  fuel  cell  is  operating  under  load,  the  data  obtained 
from  the  HP  milliohm  meter  (Model  4328A),  does  not  give  the 
values  of  the  resistance  directly,  R ,  but  only  reflects  the  change 
in  potential  drop  across  the  membrane  with  change  in  current: 


dr}  _  d(IR) 

d7  ~~  d  / 


(42) 


The  potential  drop  across  the  membrane  (i})  for  an  applied 
current  (/)  is  obtained  by  integrating  Eq.  (42).  The  area-specific 
membrane  resistance  R  (expressed  as  product  of  resistance  and 
the  area  of  the  membrane)  is  then  obtained  from  this  calculated 
potential  drop  at  each  current  density,  i: 


4.  Results  and  discussion 

4.1.  Current-voltage  plots 

We  first  discuss  the  effect  of  different  operating  parameters 
on  the  current-voltage  performance.  The  current-voltage  curves 
obtained  experimentally  are  shown  in  Fig.  3.  These  data  show  an 
overall  increase  in  performance  upon  application  of  gas  pressure 
at  either  the  anode  side  or  the  cathode  side.  The  pressure  at  the 
anode  side  slightly  enhances  the  performance.  When  pressure  is 
applied  to  the  cathode  region  the  performance  enhancement  is 


R_  >1  _  /aoiC^/d/ld/  (43) 

i  i 

4.2.2.  Validation  of  the  analytical  model  against  the 
experimental  data  and  the  numerical  solution 
4.2.2. 1.  Analytical  model  and  experimental  data.  The  mem¬ 
brane  resistance  results  (obtained  using  the  analytical  model,  Eq. 
(21)  and  the  experimental  data)  of  Nation  1 1 2  are  shown  in  Fig .  4 
as  a  function  of  the  current  density  at  70  °C.  The  parameters  used 
in  the  model  are  given  in  Table  1.  The  resistance  is  normalized 
with  the  resistance  of  the  membrane  at  open  circuit  condition 
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Fig.  4.  Comparison  of  the  analytical  solution  with  the  measured  values  for  the 
membrane  resistance  of  a  H2/O2  fuel  cell  employing  Nafion  112  membrane  and 
operated  at  70  °C  with  fully  saturated  gases  and  no  pressure  difference  between 
the  anode  and  the  cathode. 

(R  ),  for  facilitating  comparison  between  the  various  operat¬ 
ing  conditions.  In  this  study  the  resistance  R  was  found  to  be 
85  m£2  cm-2 .  The  data  shown  in  Fig.  4  corresponds  to  the  fully 
saturated  feed  condition  with  no  gas  pressure  difference  between 
both  the  anode  and  the  cathode  sides.  This  is  chosen  as  the  ref¬ 
erence  case  for  comparing  the  resistances  obtained  for  all  the 
other  operating  conditions  considered  in  this  study.  As  observed 
in  Fig.  4,  the  measured  (and  predicted)  membrane  resistance  is 
almost  constant  up  to  a  current  density  of  1.2  A  cm-2.  Beyond 


this  point  there  is  a  sharp  increase  in  the  resistance  of  the  mem¬ 
brane,  with  an  increase  in  the  current  density.  This  increase  in 
the  resistance  is  approximately  10%  (~8.5  m ST2  cm-2),  when 
the  current  density  increases  from  1.2  to  2.0  A  cm-2.  A  con¬ 
stant  value  of  the  resistance  with  has  been  reported  in  literature 

[36]  for  thin  membranes  (60  p,m).  For  Nafion  117  (~200  pun),  a 
noticeable  increase  in  the  membrane  resistance  with  an  increase 
in  the  current  density  has  been  reported  [37]  even  at  much  lower 
current  density  region  (0.2-0. 7  A  cm-2)  and  this  is  attributed  to 
the  partial  dehydration  of  the  membrane.  The  sharp  increase  in 
the  resistance  value  found  in  this  study  at  higher  current  densi¬ 
ties  indicates  the  possibility  of  partial  membrane  dehydration  at 
these  current  densities  even  for  thin  membranes  (50  |xm).  The 
analytical  model  developed  here  is  able  to  predict  this  trend  in  the 
membrane  resistance  to  a  remarkable  extent  (within  1%  error). 

The  term  in  Eq.  (15)  which  is  a  measure  of  the  electro- 
osmotic  drag,  increases  with  the  increase  in  the  current  density. 
The  difference  in  the  capillary  radius  (rc2  —  >'c  1 )  is  the  measure 
of  the  water  brought  into  the  membrane  by  the  Darcy’s  flow. 
The  greater  the  difference  in  the  capillary  radius  the  less  water 
is  brought  into  the  membrane.  The  GDL  model  (Eqs.  (36)  and 

(37) )  predicts  an  increase  in  the  difference  of  the  capillary  pres¬ 
sure  with  an  increase  in  the  current  density.  This  in  turn  results 
in  a  decrease  in  the  difference  of  the  radius  at  the  two  membrane 
boundaries.  Before  the  critical  current  density  (1.2  A  cm-2  in 
this  study)  is  reached,  the  water  loss  at  the  anode  side  of  the 
membrane  due  to  the  electro-osmotic  drag  is  compensated  by 
the  water  brought  in  by  Darcy’s  flow  resulting  in  uniform  water 
content  across  the  membrane  and  a  constant  resistance.  With  the 


Table  1 


Parameter  values  used  in  this  study 


Parameter 

Symbol 

Value/expression 

S.I.  unit 

Reference 

Cell  Temperature 

T 

353 

K 

* 

Porosity  of  the  GDL 

£ 

0.4 

- 

[40] 

Length  of  the  GDL 

Ld 

2.5  x  10-4 

m 

[41] 

Length  of  the  membrane 

Lm 

50  x  10-6 

m 

** 

Absolute  permeability  of  the  GDL 

KL 

7  x  10-6 

m2 

[40] 

Mean  radius  of  the  GDL 

rmd 

1.37  x  10~6 

m 

[40] 

Contact  angle  of  water  with  the  GDL 

eA 

2.61 

radian 

[41] 

Fraction  of  the  hydrophobic  pores  in  the  GDL 

f 

0.98 

- 

[40] 

Standard  deviation  of  the  pore  size  distribution  in  the  GDL 

<rd 

2.44  x  10~6 

m 

[40] 

Contact  angle  of  water  with  the  membrane 

em 

1.5711 

3816.44  N 

T  —  46.13  j 

radian 

[27] 

Partial  pressure  of  water 

Pw 

exp  ^1 1.6832  — 

)  x  105 

Nm"2 

[42] 

Molar  concentration  of  water 

Cw 

18 

mol  m-3 

[42] 

1.1.603  -0.00053717’ 

Viscosity  of  water 

At 

(2695.3-6.67)  > 

:  10-6 

Pas^1 

[42] 

Surface  Tension  of  water 

Y 

(0.12398-0.000173937) 

NitT1 

[41] 

Electro  osmotic  drag  coefficient  of  the  membrane 

,  55c,p  (  4000 

l  1 

PU-314 

V  303.15 

Tj) 

Faraday’s  constant 

F 

96487 

C 

[31] 

Absolute  permeability  of  the  membrane 

ism 

^abs 

5.2  x  10^20 

m2 

[43] 

Maximum  water  content  in  the  membrane 

0.44 

m3  m~3 

[33] 

Resistance  of  the  membrane  at  reference  condition 

7?a 

8.5 

m£2  cm-2 

** 

Atmospheric  pressure 

Patm 

1.023  x  105 

Nn-r2 

* 

Temperature  of  anode  humidifier 

T 

358 

K 

* 

Temperature  of  cathode  humidifier 

T 

353 

K 

Operating  conditions. 
Measured  value. 
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Fig.  5.  Comparison  of  the  analytical  solution  with  the  measured  values  for  the 
membrane  resistance  of  a  H2/O2  fuel  cell  employing  Nation  112  membrane 
and  operated  at  70  °C  with  fully  saturated  gases,  (a)  No  cathode  gas  pressure 
(analytical  model),  (b)  no  cathode  gas  pressure  (experimental  data),  (c)  cathode 
gas  pressure  of  0.5  atm  (analytical  model),  (d)  cathode  gas  pressure  of  0.5  atm 
(experimental  data),  (e)  cathode  gas  pressure  of  1 .0  atm  (analytical  model),  (f) 
cathode  gas  pressure  of  1  atm  (experimental  data),  (g)  cathode  gas  pressure  of 
2  atm  (analytical  model),  and  (h)  cathode  gas  pressure  of  2  atm  (experimental 
data). 

increase  in  the  current  density  the  electro-osmotic  drag  becomes 
predominant  thereby  causing  less  water  content  in  the  mem¬ 
brane  and  higher  resistance.  The  critical  current  density  (the 
point  of  onset  of  the  membrane  partial  dehydration),  depends 
on  the  membrane  properties  and  the  operating  conditions. 

The  resistance  of  the  membrane  as  a  function  of  current 
density  for  different  cathode  gas  pressures  is  shown  in  Fig.  5. 
These  results  show  that  there  is  an  increase  in  resistance  with  the 
increase  in  the  current  density  at  all  cathode  gas  pressure  values. 
When  the  pressure  is  2  atm  (highest  pressure  considered  in  this 
study),  there  is  a  pronounced  increase  (10%)  in  the  resistance 
with  an  increase  in  the  current  density.  But  at  the  lower  cathode 
pressure  conditions  (0.5  and  1  atm)  the  increase  in  resistance  is 
not  much  pronounced.  Another  interesting  observation  is  that 
in  all  these  three  cases  the  resistance  values  are  less  than  the 
resistance  obtained  at  the  reference  case  considered  (no  cath¬ 
ode  pressure).  All  these  observations  can  be  explained  using 
the  model  developed.  The  electro-osmotic  drag  increases  with 
increase  in  the  current  density,  causing  a  water  content  gradient 
across  the  membrane.  Darcy’s  flow  reduces  this  water  content 
gradient.  As  in  the  case  of  no  pressure  difference  across  the 
membrane,  the  increase  in  current  density  increases  the  capil¬ 
lary  pressure  gradient  across  the  membrane.  This  increase  in 
the  capillary  pressure  gradient  helps  acts  as  a  driving  force 
for  the  Darcy’s  flow.  But  at  higher  current  densities  (close  to 
2  A  cm-2)  the  electro-osmotic  drag  exceeds  the  Darcy’s  flow. 
The  application  of  a  small  pressure  at  the  cathode  acts  as  an 
additional  driving  force  and  pumps  more  water  into  the  mem¬ 
brane  at  higher  current  densities,  thus  effectively  decreasing  the 
membrane  resistance.  When  relatively  higher  pressure  (2  atm) 
is  applied  the  capillary  pressure  difference  and  hence  the  water 
content  gradient  across  the  membrane  becomes  significantly 
larger  at  higher  current  densities  and  this  is  not  compensated 


Fig.  6.  Comparison  of  the  analytical  solution  with  the  measured  values  for  the 
membrane  resistance  of  a  H2/O2  fuel  cell  employing  Nafion  112  membrane  and 
operated  at  70  °C  with  fully  saturated  gases,  (a)  No  anode  gas  pressure  (analytical 
model),  (b)  no  anode  gas  pressure  (experimental  data),  (c)  anode  gas  pressure 
of  0.5  atm  (analytical  model),  (d)  anode  gas  pressure  of  0.5  atm  (experimental 
data),  (e)  anode  gas  pressure  of  2.0  atm  (analytical  model),  and  (f)  anode  gas 
pressure  of  2  atm  (experimental  data). 

by  the  additional  convective  driving  force  provided  by  the  dif¬ 
ference  in  gas  pressure. 

The  effect  of  anode  gas  pressure  on  the  membrane  resis¬ 
tance  at  70  °C  is  depicted  in  Fig.  6.  We  observe  a  decrease  in 
the  membrane  resistance  (approximately  17%)  with  an  increase 
in  the  current  density  from  0.1  to  2  A  cm-2.  This  decrease  in 
the  membrane  resistance  with  increase  in  the  current  density 
is  nearly  the  same  for  both  low  and  high  anode  gas  pressures 
considered  in  this  study.  When  we  analyze  the  situation  using 
the  model,  we  find  that  with  the  application  of  pressure  at  the 
anode  the  capillary  pressure  at  the  anode  side  of  the  membrane 
increases  resulting  in  increase  in  water  content  on  the  anode  side. 
This  complements  the  Darcy’s  flow  in  compensating  the  water 
loss  due  to  the  electro-osmotic  drag,  ensuring  low  membrane 
resistance. 

Fig.  7  shows  the  resistance  of  the  membrane  as  a  function 
of  current  density  when  dry  O2  gas  is  used  as  the  oxidant.  In 
this  case,  the  resistance  decreases  gradually  from  about  36% 
more  than  the  corresponding  reference  case  value  to  almost  equal 
value  to  that  of  the  reference  case  when  the  current  density  is 
increased  from  0.1  to  2  A  cm-2.  In  the  lower  current  density 
region  the  water  flux  by  Darcy’s  flow  is  comparatively  less  due  to 
less  water  content  gradient  across  the  membrane.  Thus  the  water 
dragged  into  the  PEM  by  the  proton,  which  only  depends  on  the 
operating  current  density  is  higher  than  the  flow  of  water  from  the 
cathode  to  the  anode,  giving  rise  to  a  higher  value  of  resistance 
than  the  reference  case.  This  loss  is  gradually  compensated  by 
the  increased  production  of  water  at  the  cathode  which  increases 
the  water  content  gradient  across  the  membrane  aiding  in  more 
flux  of  water  into  the  membrane  due  to  the  Darcy’s  flow. 

Fig.  8  shows  the  membrane  resistance  for  the  case  of  dry 
H2  gas  and  fully  saturated  cathode  gas  with  no  pressure  differ¬ 
ence  between  the  anode  and  the  cathode.  At  all  current  densities 
the  resistance  is  larger  than  the  reference  case.  This  can  be 
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Fig.  7.  Comparison  of  the  analytical  solution  with  the  measured  values  for  the 
membrane  resistance  of  a  H2/O2  fuel  cell  employing  Nafion  112  membrane 
and  operated  at  70  °C  with  no  pressure  difference  between  the  anode  and  the 
cathode,  (a)  Fully  humidified  gases  (analytical  model),  (b)  fully  humidified  gases 
(experimental  data),  (c)  dry  cathode  gas  (O2)  (analytical  model),  and  (d)  dry 
cathode  gas  (O2)  (experimental  data). 

explained  by  the  unequal  compensation  of  water  loss  due  to 
electro-osmosis  by  the  Darcy’s  flow.  Water  is  brought  into  the 
membrane  from  the  cathode  by  Darcy’s  flow  and  the  proton 
drags  water  from  the  anode  side  along  with  it  toward  the  cath¬ 
ode  (electro-osmotic  drag).  Since  there  is  no  water  present  at  the 
anode  membrane  boundary,  the  water  brought  in  by  the  Darcy’s 
flow  is  dragged  back  towards  the  cathode,  resulting  in  severe 
dehydration  of  the  anode  side  of  the  membrane  and  hence  high 
membrane  resistance.  This  effect  is  significantly  reduced  with 
an  increase  in  the  current  density  as  more  and  more  water  is 
produced  at  the  cathode  and  relatively  more  water  is  dragged 
by  Darcy’s  flow  from  the  cathode  side  to  the  anode.  This  can 
be  observed  from  Eq.  (21).  With  increase  in  current  density  the 
capillary  radius  decreases  sharply  and  thereby  increasing  the  dif- 


Fig.  8.  Comparison  of  the  analytical  solution  with  the  measured  values  for  the 
membrane  resistance  of  a  H2/O2  fuel  cell  employing  Nafion  112  membrane 
and  operated  at  70  °C  with  no  pressure  difference  between  the  anode  and  the 
cathode,  (a)  Fully  humidified  gases  (analytical  model),  (b)  fully  humidified  gases 
(experimental  data),  (c)  dry  anode  gas  (H2)  (analytical  model),  and  (d)  dry  anode 
gas  (H2)  (experimental  data). 


Fig.  9.  Comparison  of  numerical  solutions  with  the  measured  values  and  the 
analytical  solution  for  the  membrane  resistance  of  a  H2/O2  fuel  cell  employing 
Nafion  112  membrane  and  operated  at  70  °C  with  fully  humidified  gases,  (a)  No 
gas  pressure  at  the  anode  and  the  cathode  (analytical  model),  (b)  no  gas  pressure 
at  the  anode  and  the  cathode  (experimental  data),  (c)  no  gas  pressure  at  the  anode 
and  the  cathode  (numerical  model),  (d)  cathode  gas  pressure  of  2  atm  (analytical 
model),  (e)  cathode  gas  pressure  of  2  atm  (experimental  data),  (f)  cathode  gas 
pressure  of  2  atm  (numerical  model),  (g)  anode  gas  pressure  of  2  atm  (analytical 
model),  (h)  anode  gas  pressure  of  2  atm  (experimental  data),  and  (i)  anode  gas 
pressure  of  2  atm  (numerical  model). 

ference  in  the  capillary  radius  between  anode  and  cathode  side, 
causing  the  resistance  to  drop. 

Hence  the  optimum  condition  for  fuel  cell  from  the  point  of 
view  of  membrane  performance,  is  operating  the  system  with  no 
or  minimum  humidity  of  the  cathode  gas,  with  the  application 
of  a  small  pressure  at  the  cathode  (around  0.5  atm)  and  com¬ 
paratively  large  pressure  (around  2  atm)  at  the  anode  with  fully 
humidified  anode  gas. 

4. 2. 2.2.  Comparison  of  analytical,  numerical  solution  with 
experimental  data.  Fig.  9  presents  a  comparison  of  the  numer¬ 
ical  solution  (Eq.  (18)  along  with  Eq.  (23))  with  the  analytical 
solution  (21)  and  the  experimental  data  for  three  different  oper¬ 
ating  conditions.  It  is  shown  that  when  there  is  no  pressure 
difference  between  the  anode  and  the  cathode  the  numerical 
solution  and  the  analytical  approximation  are  identical.  But 
when  the  pressure  difference  exists,  the  numerical  solution  is 
relatively  close  (error  is  less  than  1%)  to  the  experimental  data, 
while  the  analytical  solution  differs  by  3%  from  the  experimental 
data.  But  this  small  difference  should  not  prevent  us  from  taking 
advantage  of  the  analytical  solution.  The  approximate  analytical 
solution  can  reduce  the  computational  time  in  a  full-scale  model. 

5.  Conclusion 

The  resistance  of  Nafion  112  membrane  was  investigated 
under  different  operating  conditions.  It  was  observed  that  though 
application  of  pressure  at  both  the  anode  side  and  the  cathode 
side  helps  in  increasing  the  performance  of  the  cell,  pressure  at 
the  anode  helps  to  reduce  the  membrane  resistance  to  a  larger 
extent.  The  effect  of  dry  and  fully  humidified  gases  on  the 
membrane  performance  was  studied.  It  was  observed  that  anode 
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humidity  decreases  the  membrane  performance.  The  optimum 
operating  pressure  conditions  were  predicted  using  the  analyt¬ 
ical  expression  developed  for  the  membrane  resistance.  Also  a 
simple  model  to  predict  the  membrane  resistance  was  devel¬ 
oped  which  was  able  to  predict  the  trend  in  the  PEM  resistance 
with  good  accuracy  and  also  provides  certain  physical  insight 
for  better  membrane  performance. 


with  water  while  in  hydrophobic  pores  the  pores  above  the  criti¬ 
cal  radius  are  filled  with  water.  This  concept  lets  us  fix  the  limits 
of  integration  in  both  the  case.  The  approach  followed  here  is 
similar  to  that  of  White  et  al.  [39].  Here  the  pore  size  distribution 
function  is  assumed  to  be  uniform  for  both  the  type  of  pores. 

B.l.  Derivation  of  water  content  of  the  GDL 
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Appendix  A 

A.l.  Calculation  of  water  flux  across  the  membrane 


The  flux  of  water  in  the  anode  (same  as  the  flux  of  water  at 
the  inlet  of  anode  GDL)  is  given  by 


Nm/d  -  — 
lyW  —  -  - 


IF  \  1 


(A.l) 


where  u>a  is  the  moles  of  water  leaving  the  anode  humidifier.  wa 
at  the  operating  temperature  of  the  anode  humidifier  (75  °C)  is 
numerically  equal  to  0.3743.  Hence  at  the  flux  of  water  in  the 
GDL  is 


i !f\ir)  —  f  a(r)dr  (B.l) 

Jrc 

where  the  pore  size  distribution  function  is  assumed  to  be  (based 
on  the  fit  obtained  from  typical  pore  size  distribution  data  for 
the  GDL  material  [40]): 


1 

a(r)  —  - ^=exp 

ra^/2n 


In  r  —  In  rnr 


V2c 


(B.2) 


Hence  upon  integration: 
f\(r'c)  — 


erf^c) 


(B.3) 


where  r’c  is  ((lnrc  —  lnrmd) / \fla) 

Pore  size  distribution  for  hydrophobic  gas  pores  is  given  by 


Nm’a 

tvw 


0.31 

~T~ 


(A.2) 


The  flux  of  water  due  to  electro-osmotic  drag  across  the  mem¬ 
brane  is  given  by 


drag 


W 


F 


(A. 3) 


For  the  Nafion  membrane  equilibriated  with  liquid  water  the 
electro-osmotic  coefficient  is  given  as  2.55  [27].  Hence  the  flux 
of  water  due  to  electro-osmotic  drag  is 


..drag  2.55  i 

Nw  6  ~  - 

W 


(A.4) 


At  all  operating  current  densities,  the  flux  due  to  the  electro- 
osmotic  drag  (A.4)  is  greater  than  the  flux  of  water  at  the  anode 
side  of  the  membrane  (A.2).  This  ensures  that  the  water  brought 
into  the  membrane  from  the  cathode  does  not  flow  into  the  anode, 
but  only  humidifies  the  anode  side  of  the  membrane.  Thus  the 
net  water  flux  across  the  membrane  is  negligible  at  the  current 
operating  conditions. 


flir)  —  [  a(r)dr  (B.4) 

Jo 

where  the  pore  size  distribution  function  is  assumed  to  be  the 
same  log  normal  distribution. 

Now  we  get: 


020c) 


1  +  erf (Fc) 
2 


(B.5) 


The  fraction  of  the  hydrophilic  pores  filled  with  gas  or  the 
effective  hydrophilic  porosity  is 


01  =  s(l  -  f)Mr'c) 


(B.6) 


where  e  is  the  porosity  of  the  GDL  and  /  is  the  fraction  of 
hydrophobic  pores  in  the  GDL  media.  And  for  hydrophobic 
pores: 

02  =  efihtfc)  (B.7) 

Total  fraction  of  pores  filled  with  gas  is 


Appendix  B 


0T  =  01  +  02  =  1 


(B.8) 


The  GDL  is  considered  to  be  a  porous  medium  with  both 
hydrophilic  and  hydrophobic  pores.  The  relationship  between 
the  capillary  radius  and  water  content  is  derived  based  on  the 
concept  of  wetting  [38].  The  critical  radius  required  for  water 
to  wet  the  pores  depends  on  the  nature  of  the  pore.  In  the  case 
of  hydrophilic  pores,  the  radii  below  the  critical  radius  are  filled 


Hence  the  total  gas  filled  pores,  is  given  by 


0T  = 


1  -  erf(^)(l  -  202) 
2 


(B.9) 


The  liquid  water  saturation  or  the  water  content  of  the  medium 
is  the  fraction  of  pores  filled  with  water  which  is  1  —  0t  and 
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substituting  for  rc  we  have: 

Sw  —  l  ^1  +  0  -  202)erf 

l ln(— ((2y  cos  9)/ pc))  -  ln(rmd) 

x  { 

This  expression  is  used  in  the  GDL  model  for  relating  the  GDL 
properties  and  the  water  content  in  the  medium. 
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